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^ ! ABSTRACT 

a • 

1 I ' Using Planck data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions 
, by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at 
, the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow 

^ ' amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to 
' the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to 72 ± 60kms"'. This constitutes less than 
, 1 % of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. From a subset of this cluster sample Planck finds 
, the radial peculiar velocity rms to be below 800km s"' at the 95 % confidence level, which is around three times the ACDM prediction for the 
' typical cluster radial velocity rms at z = 0. 15. Planck data also set strong constraints on the local bulk flow in volumes centred on the Local Group. 
' There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind 
, search for bulk flows in this sample has an upper limit of 254 km s"' (95 % confidence level) dominated by CMB confusion and instrumental 
noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, Planck is 
I ■ able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The Planck constraints on peculiar 

. • ' velocities and bulk flows are thus consistent with the ACDM scenario. 

> ■ 

• i-H ' Key words, cosmology: observations - cosmic microwave background - large-scale structure of the Universe - galaxies: clusters: general 

^ : 

- ■ ■ 1 . Introduction of dark matter kept growing gravitationally. It was only after the 

epoch of hydrogen recombination that the baryons, which had 

Today we have a cosmological model that appears to fit all just decoupled from the CMB, could freely fall into the potential 

available data. Nevertheless, it is important to continue to test ^^jj^ j^^,^ ^^^^^^ component, 
this picture. Peculiar velocities provide an important way to do 

so. According to the standard ACDM scenario, gravity drives Since then, the gravitational infall of matter into potential 

the growth of inhomogeneities in the matter distribution of the wells has been conditioned by the density field and the uni- 

Universe. After the radiation-matter equality epoch, fluctuations versal expansion rate. On large scales, where baryonic physics 

in the dark matter component were largely unaffected by the and non-linear evolution may safely be neglected, the continuity 

Thomson interaction binding the evolution of baryons and pho- equation provides a simple link between the matter density field 

tons of the cosmic microwave background radiation (CMB). and the peculiar velocity field. In particular, in a ACDM sce- 

During that epoch, the inhomogeneities in the spatial distribution nario, it predicts that peculiar velocities must show typical cor- 

relation lengths of the order of 20-40 /i ' Mpc, and their growth 

* Corresponding author: C.Hemandez-Monteagudo: chm@cefca.es must have practically frozen since the onset of the accelerated 
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expansion (see the review of, e.g., Strauss & Willick 1995). By 
averaging the peculiar velocity field on scales corresponding to 
galaxy groups and clusters today it is possible to obtain linear 
theory predictions for the root mean square (rms) of the radial 
peculiar velocity of those structures. These predictions typically 
amount to about 230 km s"', (see, e.g., Hemandez-Monteagudo 
& Sunyaev 2010), with a weak dependence on the galaxy clus- 
ter/group mass. If instead one looks at the velocity amplitude 
for extended or correlated motion of matter on larger scales, one 
finds that it decreases when larger volumes are considered, but 
should still be at the level of 50-100 kms ' for radii of a few 
hundred Mpc (see, e.g., figure 2 in Mak et al. 2011). The de- 
tection of these large-scale, coherent flows of matter (hereafter 
referred to as bulk flows) has been the subject of active investi- 
gation for several decades (e.g., Tonry & Davis 1981; Aaronson 
et al. 1982; Dressier et al. 1987; Dekel et al. 1993; Lauer & 
Postman 1994; Hudson et al. 1999; Willick 1999; Riess 2000, 
to cite just a few). One crucial problem that most of those works 
encounter is related to the need to accurately determine distances 
to galaxies in order to subtract the Hubble flow-induced velocity. 

During the nineties Lauer & Postman ( 1994), Wiflick (1999), 
and Hudson et al. (1999) claimed the existence of large scale 
bulk flows with amplitudes of 350-700 km s"' in local spheres of 
radii 60-150 /i"' Mpc, with somewhat discrepant directions. At 
the turn of the millennium, however, Riess (2000) and Courteau 
et al. (2000) reported the lack of any significant local bulk flow 
up to depths of about 150/2 'Mpc, in apparent contradiction 
with the previous works. More recently, claims of the presence of 
a large-scale, large amplitude peculiar velocity dipole have been 
raised again by various authors. While some works (Hudson 
et al. 2004; Watkins et al. 2009; Feldman et al. 2010) find ev- 
idence for a peculiar local velocity dipole of about 400 kms"' 
(and in tension with ACDM predictions), others find lower am- 
plitudes for the local bulk flow, (e.g., Nusser & Davis 201 1; Ma 
& Scott 2012). 

At much larger depths (up to z ~ 0.2-0.3) there are also 
claims (Kashlinsky et al. 2008, 2010; Abate & Feldman 2012) 
of yet higher amplitude bulk flows (~ 1 000-4000 km s"' ). These 
cannot be accommodated within a ACDM context, since the the- 
ory predicts bulk flows to be negligible on the very largest scales. 
Moreover, these results are in contradiction with other works, 
(e.g., Keisler 2009; Osborne et al. 201 1 ; Mody & Hajian 2012). 

Some of the most recent results on bulk flows extending to 
large distances are based on the study of the kinetic Sunyaev- 
Zeldovich effect (hereafter kSZ; Sunyaev & Zeldovich 1980). 
This effect is due to the Doppler kick that CMB photons ex- 
perience in Thomson scattering off free electrons moving with 
respect to the CMB rest frame. This process introduces inten- 
sity and polarization anisotropics in the CMB along the direc- 
tion of massive clouds of ionized material, such as galaxy clus- 
ters and groups, but produces no distortion of the CMB spec- 
trum. The kSZ effect has been theoretically exploited to char- 
acterize the growth of velocity perturbations (e.g.. Ma & Fry 
2002; Hemandez-Monteagudo et al. 2006; Zhang et al. 2008), 
to search for missing baryons (DeDeo et al. 2005; Hemandez- 
Monteagudo & Sunyaev 2008; Ho et al. 2009; Hemandez- 
Monteagudo & Ho 2009; Shao et al. 2011), and to study bulk 
flows in the local Universe (Kashlinsky et al. 2008, 2010; Keisler 
2009; Osborne et al. 2011; Mak et al. 2011; Mody & Hajian 
2012; Lavaux et al. 2012). Very recently. Hand et al. (2012) 
claimed a detection of the kSZ effect when combining spec- 
troscopic galaxy data from the Baryonic Acoustic Oscilation 
Survey (BOSS) with CMB data from the Atacama Cosmology 
Telescope (ACT), after searching for the kSZ pairwise mo- 



mentum (e.g., Groth et al. 1989; Juszkiewicz et al. 1998). On 
sub-cluster scales, as predicted by, e.g., Inogamov & Sunyaev 
(2003), some weak evidence for kSZ has also been reported by 
the MUSTANG experiment (Mroczkowski et al. 2012). 

In this paper we shall focus on the constraints that Planck^ 
can set on the kSZ-induced temperature anisotropics. These are 
given by the line-of-sight integral 

ST f Ve ■ n 

—(«) = - dla-jHe , (1) 

To J c 

where ctj is the Thomson scattering cross-section, is the phys- 
ical electron number density, Ve denotes the electron peculiar ve- 
locity, c the speed of light, and h the direction of observation on 
the sky. We are adopting here a reference frame centred on the 
observer's position, and hence infalling electrons will have nega- 
tive radial velocities. Note that, unlike in other approaches based 
upon galaxy redshift surveys, the distance to the cluster is irrel- 
evant to its peculiar velocity estimation. Since the expected kSZ 
signal coming from an individual cluster is smaller than the typi- 
cal level of intrinsic CMB temperature fluctuations, we shall ap- 
ply various filters which attempt to minimize the impact of other 
signals on the angular positions of a sample of galaxy clusters, 
and use these to extract statistical constraints on the kSZ signal 
in those sources. In the standard ACDM scenario one expects to 
have matter at rest with respect to the CMB on the largest scales, 
and hence roughly the same number of clusters with positive and 
negative radial velocities. This means that the mean or monopole 
of kSZ estimates should be consistent with zero, although there 
are inhomogeneous scenarios (addressed in Sect. 4.3) in which 
the average velocity of clusters may differ from zero. Likewise it 
is possible to set constraints on the kSZ-induced variance in the 
CMB temperature anisotropics measured along the direction of 
galaxy clusters. This is a direct probe of the rms peculiar velocity 
of those objects with respect to the CMB, and can be compared 
to theoretical predictions. In this context, it has been mentioned 
above that the motion of matter is predicted to occur in bulk 
flows with coherence on scales of about 30 /i 'Mpc. If these bulk 
flows are local and the observer is placed inside them, then they 
should give rise to a dipolar pattern in the kSZ measurements of 
individual clusters (Kashlinsky & Atrio-Barandela 2000). If they 
are instead distant, then the projection of the coherence length of 
kSZ measurements on the sky should shrink down to a few de- 
grees (Hemandez-Monteagudo et al. 2006). Therefore it is pos- 
sible to use the set of individual kSZ estimates from galaxy clus- 
ters to place constraints on the monopole (mean), variance, and 
dipole of the peculiar velocities of the cluster population. While 
some of our statistical tools target the kSZ signal in each clus- 
ter separately, others are particularly designed to probe the local 
bulk blow and set constraints on the kSZ dipole at the positions 
of clusters, as will be shown below. 

This paper is organized as follows. In Sect. 2 we describe 
the data used, both for the CMB and large-scale structure. The 
statistical tools we use for the kSZ detection are described in 
Sect. 3, and the results obtained from them are presented in 
Sect. 4. We examine the robustness of our results in Sect. 5. 
Finally, in Sect. 6, we discuss the cosmological implications 
of our findings and conclude. Throughout this paper, we use a 



' Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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cosmological parameter set compatible with WMAP-1 observa- 
tions (Komatsu et al. 201 1): density parameters = 0.272 and 
- 0.728; Hubble parameter h - 0.704; 8 /r'Mpc normaliza- 
tion (Tg = 0.809; and «s — 0.96 for the spectral index of scalar 
perturbations. 

2. Data and simulations 

2.1. Planck data 

Planck (Tauber et al. 2010; Planck Collaboration I 2011) is the 
third generation space mission to measure the anisotropy of the 
CMB. It observes the sky in nine frequency bands covering 30- 
857 GHz with high sensitivity and angular resolution from 3 1 ' 
to 4.39'. The Low Frequency Instrument (LFI; Mandolesi et al. 
2010; BersanelH et al. 2010; Mennella et al. 2011) covers the 
30, 44, and 70 GHz bands with amplifiers cooled to 20 K. The 
High Frequency Instrument (HFI; Lamarre et al. 2010; Planck 
HFI Core Team 2011a) covers the 100, 143, 217, 353, 545, and 
857 GHz bands with bolometers cooled to 0. 1 K. Polarization is 
measured in all but the highest two bands (Leahy et al. 2010; 
Rosset et al. 2010). A combination of radiative cooling and 
three mechanical coolers produces the temperatures needed for 
the detectors and optics (Planck Collaboration II 2011). Two 
data processing centres (DPCs) check and calibrate the data and 
make maps of the sky (Planck HFI Core Team 2011b; Zacchei 
et al. 2011). Planck's sensitivity, angular resolution, and fre- 
quency coverage make it a powerful instrument for Galactic 
and extragalactic astrophysics as well as cosmology. Early as- 
trophysics results are given in Planck Collaboration VIII-XXVI 
201 1, based on data taken between 13 August 2009 and 7 June 
2010. Intermediate astrophysics results are now being presented 
in a series of papers based on data taken between 13 August 
2009 and 27 November 2010. 

Although the 70 GHz LFI channel was included initially 
in the analysis, it was found that constraints were practically 
identical when using HFI frequency maps alone (see details in 
Table 1). Measuring the kSZ effect requires avoiding, in the 
best possible way, contamination by the much stronger thermal 
Sunyaev-Zeldovich effect (hereafter tSZ; Sunyaev & Zeldovich 
1972). While in theory observations at 217 GHz, close to the 
zero of the tSZ emission, should not suffer much from tSZ 
contamination, it is necessary to account for the broad spec- 
tral band of each detector and each channel, (Planck HFI Core 
Team 201 lb). In terms of the tSZ effect, the 'effective' frequen- 
cies of the HFI channels (i.e., those frequencies at which the 
tSZ emission is equal to its integral over the frequency band) 
are listed in the second column of Table 1 . Raw HFI frequency 
maps are useful for testing for systematic effects associated with 
foreground emission, tSZ spectral leakage, or FWHM charac- 
terization errors. Note also that HFI frequency maps are pro- 
duced in thermodynamic temperature units, so that both pri- 
mary CMB and kSZ emission have constant amplitude across 
frequency channels. In the third column we display the corre- 
sponding tSZ Comptonization parameter (ysz) to hJ conversion 
factors. The Comptonization parameter ysz is a dimensionless 
line-of-sight integral of the gas pressure, 

ysz = 1 dlo-Tn^'^^^, (2) 

J nieC" 

with Tf. and nif^ the electron temperature and rest mass, and k^ 
the Boltzmann constant, (Sunyaev & Zeldovich 1972). 

Two different strategies are used in this paper to measure 
the kSZ effect at the positions of the cluster catalogue. The first 



Table 1. Nominal and tSZ-effective frequencies, Ar to ysz con- 
version factors and FWHMs for each HFI channel used in this 
paper. The second column provides the effective frequency of 
the channels after considering the spectral dependence of the 
non-relativistic tSZ effect and the finite response of the HFI de- 
tectors. The third column displays the conversion factor between 
SZ Comptonization parameter (ysz) and thermodynamic temper- 
ature (in K). The fourth column provides, for each HFI channel, 
the average FWHM value of the effective Gaussian beam at map 
level, as described in Planck HFI Core Team (201 lb). 



HFI nominal HFI effective vsz/Ar FWHM 



frequency 


frequency 






[GHz] 


[GHz] 




[arcmin] 


100 


103.1 


-0.2481 


9.88 


143 


144.5 


-0.3592 


7.18 


217 


222.1 


5.2602 


4.87 


353 


355.2 


0.1611 


4.65 


545 


528.5 


0.0692 


4.72 


857 


775.9 


0.0380 


4.39 



consists of estimating directly the kSZ signal at MCXC clus- 
ter positions from the original HFI frequency maps, using both 
aperture photometry and matched multi-band filtering. The sec- 
ond consists of first producing a map of the CMB and kSZ effect 
that is nearly free from tSZ contamination before estimating the 
kSZ emission from MCXC clusters using the aperture photome- 
try and single frequency matched filtering. As described below, 
this map makes use of both HFI and LFI data. 

2.1 .1 . The two-dimensional internal linear combination map 

In the absence of a fully reliable model of foreground emission 
(including number of foregrounds, emission laws, and coherence 
of their emission across Planck frequencies), a minimum vari- 
ance map of CMB emission can be obtained by the so-called 
"Internal Linear Combination" approach (hereafter ILC). The 
CMB map is obtained from a linear combination of input ob- 
servations, subject to the constraint that the CMB is preserved. 
I.e., for CMB -calibrated maps (in thermodynamic units) x,, the 
CMB is obtained as 2, with Yji^i - 1' the latter condition 
guaranteeing the preservation of the CMB signal. This obviously 
also preserves the kSZ signal, which has the same frequency de- 
pendence. 

The minimum variance map, however, is not necessarily that 
of minimum contamination by any particular foreground. In our 
present analysis, the measured map of CMB-nkSZ will be further 
processed, first being filtered on the basis of predicted kSZ clus- 
ter shapes and locations, to suppress contamination by the larger 
scale primary CMB, and then stacked to combine the measure- 
ments of all individual clusters. While this filtering and stacking 
will reduce contamination of the measurement by independent 
foregrounds such as Galactic dust emission, as well as by de- 
tector noise, tSZ residuals are likely to add-up coherently and 
contaminate the measurement significantly. The ILC must then 
be modified to ensure that instead of the total variance of the map 
being minimized, the contamination by tSZ must be minimized 
instead. 

It is possible to extend the ILC method to add a constraint 
to reject the tSZ effect specifically, and thus make sure that 
the CMB+kSZ map is completely free from tSZ contamination. 
The idea is similar to that used in the unbiased Multifrequency 
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Matched Filter approach (Herranz et al. 2005; Mak et al. 2011), 
i.e., a constraint is added to null the tSZ contribution to the out- 
put map. The method is described and validated on realistic sim- 
ulations in Remazeilles et al. (2011a). 

Note that this method is a special case of a multi-dimensional 
generalization of the ILC (Remazeilles et al. 2011b), in which 
several components of interest with known emission laws, can 
be recovered simultaneously with vanishing contamination from 
each other Here we consider two specific components, one with 
the CMB emission law, which comprises both primary CMB and 
kSZ, and one with the emission law of the tSZ effect (neglecting 
relativistic corrections). We refer to the map obtained by this 
method as a two-dimensional ILC (hereafter 2D-ILC). 

In detail, the 2D-ILC map used in the present analysis is ob- 
tained from all LFI and HFI maps as follows. For each frequency 
band, point sources detected by a Mexican Hat Wavelet filter at 
more than 5 cr at that frequency are masked. The masked region 
has a radius of three times the standard deviation of the Gaussian 
beam (i.e., 1.27xFWHM). The masked regions are filled in by 
interpolation using neighbouring pixels. Maps are then anal- 
ysed on a frame of spherical needlets for implementation of the 
ILC in needlet space, in a very similar way to what has been 
done on WMAP data by Delabrouille et al. (2009). However, the 
covariance matrices associated with the filter, instead of being 
computed using average covariances of needlet coefficients over 
HEALPix- superpixels, are computed from products of maps 
of needlet coefficients, smoothed using a large Gaussian beam, 
similarly to what was done by Basak & Delabrouille (2012) on 
WMAP 7-year data. The constrained ILC filter implemented is 
that of Eq. 20 of Remazeilles et al. (201 la). Thus, the exact lin- 
ear combination used to reconstruct the CMB +kSZ map depends 
both on the sky region and on the angular scale. In particular, 
on scales smaller than some of the Planck LFI and HFI beams, 
the relative weights of the corresponding lower frequency chan- 
nels become negligibly small, due to their low resolution. The 
final CMB-HkSZ map is reconstructed at 5' resolution. In order 
to carry this out, at the very smallest scales the CMB-i-kSZ map 
is reconstructed mostly from observations in the frequency chan- 
nels at 217 GHz and above. At intermediate scales (of order 10'), 
however, measurements from all HFI channels are used to recon- 
struct the final map. 

The ILC (classical or 2D version) assumes the emission law 
of the component of interest to be known. This knowledge is nec- 
essary to ensure the preservation of the signal of interest (here, 
the kSZ effect) and, for the 2D-ILC, to reject the contaminat- 
ing signal (here, the tSZ). As discussed in Dick et al. (2010), 
imperfect knowledge of the emission law can result in a signifi- 
cant loss of CMB power In practice, the effective emission laws, 
as observed by the detectors, depend on the calibration of the 
observations in each frequency channel. For the CMB and kSZ 
signals, it hence depends on the absolute calibration of all the 
Planck channels used in the analysis (here, the HFI channels). 
For the tSZ, it also depends on the accuracy of the knowledge 
of the spectral bands, and on the validity of the non-relativistic 
approximation for tSZ emission. 

For Planck HFI channels, the absolute calibration error is es- 
timated (conservatively) to be less than about 0. 1 % in the chan- 
nels calibrated on the CMB itself with the CMB dipole, and on 
CMB anisotropics themselves for cross-calibration, and less than 
a few percent on channels calibrated on the dust emission mea- 
sured by FIRAS (545 and 857 GHz channels, Planck HFI Core 
Team 2011b). Small uncertainties on the frequency dependence 



^ http://healpix.jpl.nasa.gov 



of the (CMB+ kSZ) signal may induce a large bias in the cali- 
bration of the output of the ILC. This eff'ect can be strong in the 
high SNR regime (Dick et al. 2010), which is the case here be- 
cause the strong CMB anisotropy signal itself contributes to the 
total signal. With CMB calibration errors of 0.1 %, we check on 
simulated data sets generated with the the Planck Sky Model 
(Delabrouille et al. 2012) that the corresponding error on the 
final map is small (less than 1 %). This is also confirmed on 
the actual Planck data by comparing the power spectrum of the 
CMB+kSZ map with the current CMB best fit Cf, since any se- 
rious loss of power would be immediately visible in the power 
spectrum of the reconstructed CMB+kSZ map. 

Errors in the assumed tSZ emission law (by reason of rel- 
ativistic corrections, colour correction, or mis-calibration), can 
also potentially result in residual contamination by tSZ in the 
CMB+kSZ ILC map. Note, however, that the two-dimensional 
ILC does not amplify the contamination by a mis-calibrated tSZ 
component. Uncertainties of a few percent on the tSZ frequency 
dependence (the typical size of relativistic corrections to the 
thermal SZ effect) will hence not impact the reconstruction of 
the CMB+kSZ signal by more than a few percent of the original 
tSZ. 

2.2. Tracers of moving baryons: X-ray MCXC clusters 

We use the Meta Catalogue of X-ray detected Clusters of galax- 
ies (MCXC), a compilation of all publicly available ROSAT 
All Sky Survey-based samples (NORAS, Bohringer et al. 2000; 
REFLEX, Bohringer et al. 2004; BCS, EbeUng et al. 1998, 
Ebeling et al. 2000; SGP, Cruddace et al. 2002; NEP, Henry 
et al. 2006; MACS, Ebeling et al. 2007; and CIZA, Ebeling 
et al. 2010, Kocevski et al. 2007), and serendipitous cluster cat- 
alogues (160SD, Mullis et al. 2003; 400SD, Burenin et al. 2007; 
SHARC, Romer et al. 2000, Burke et al. 2003; WARPS, Perlman 
et al. 2002, Horner et al. 2008; and EMSS, Gioia & Luppino 
1994 Henry 2004). The information was systematically homog- 
enized and duplicate entries were carefully handled, yielding a 
large catalogue of approximately 1750 clusters. The MCXC is 
presented in detail in Piffaretti et al. (2011), and has been used 
in previous Planck studies, (e.g., Planck Collaboration X 201 1). 

For each cluster the MCXC provides, among other quanti- 
ties, the position, redshift, and mass of each cluster The masses 
are estimated from the luminosities thanks to the REXCESS 
Lgoo-Msoo relation of Pratt et al. (2009). Hereafter, all cluster 
quantities with the subscript "500" are evaluated at the radius 
(/?50o) at which the average density equals 500 times the critical 
density at the cluster's redshift. In this way, Msoo is defined as 
M5Q0 = (47r/3)500pc(z)^5oo' where Pc(z) is the critical density 
at the cluster redshift z. 

For the measurement of velocities, we also need cluster opti- 
cal depths. Our approach here is based upon the study of Amaud 
et al. (2010): using REXCESS data, we either use the universal 
pressure profile and then divide by the average temperature pro- 
file to estimate a density profile, or fit directly an average density 
profile. For this purpose, we make use of Y{x), the volume inte- 
gral of gas pressure up to a radius given by x = r/R^oo: 

5'sph(-x) = dr o-jneir ) (3) 

Jo '"eC^ 

Ycy\(x) = dx I dzo-jn^r) —. (4) 

Jo J-00 '«eC 

These equations describe the spherical and cylindrical integrals 
of pressure, respectively. The vector centred on the cluster r' is 
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decomposed into a vector on the plane of the sky, x' , and a vector 
normal to this plane (given by the component z'), r' - (x' ,7l). 
With this, we use the Y(x) vs. M500 relations in Amaud et al. 
(2010), i.e.. 



}'sph(x) = F500/(X), 
FcyK-^) = i'500/(x), 

in which 

7500 = 1-38 x 10-^ £2/3(^)|- 



(5) 
(6) 



^5500 / 



X 



SOO/i^o'Mpc 



-2 



/17Q arcmin", 



(7) 



and /(x) and 7(x) are functions expressing the spheri- 
cal/cylindrical integrals of pressure around the cluster's centre, 
respectively. The factor /ivo denotes the Hubble reduced param- 
eter in units of 70kms"' Mpc"'. An observer's angular aperture 
Q translates into an effective cluster radius r — 9Da{z), with 
Da(z) the angular diameter distance to redshift z- Note that, as 
expressed above, the cylindrical case considers a full integration 
along the line of sight up to a given aperture on the plane of the 
sky, as is the case for real observations. The spherical case, in- 
stead, integrates out to a given radius in all directions, and differs 
from the cylindrical case in a geometric factor The constants in 
Eq. 7 are B500 = 3 x lO"* /i^ijMo and ay = 1.78, while E{z) is 
the Hubble parameter normalized to its current value. 



(8) 



The functions I(x) and J{x) depend on the particular model 
adopted for the pressure profile, which, in our work, is taken 
to follow the universal scaling provided in Eqs. 11 and 12 of 
Arnaud et al. (2010), for which I(x = 1) = 0.6541 and J{x = 
1) = 0.7398. As will be addressed in Sect. 5, results do not 
change significantly when adopting different choices for Arnaud 
et al. (2010) type pressure profiles, but their uncertainty is dom- 
inated by our ignorance of the gas density profile in the clusters' 
outskirts. 

In the isothermal case, the clusters' optical depth integrated 
up to xR^Qo is equal to the corresponding Y(x), modulo a 
kT linisC") factor, where f is the average spectroscopic temper- 
ature measured in a fraction of the volume enclosed by /?5oo. To 
account for this, we use the f-Msoo relation given in Arnaud 
et al. (2005): 



£(z)M5, 
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(9) 



with ax - 1-49 and A500 - 4.10 x lO'^^/z^g Mq. In this simple 
case, it is possible to derive an expression for the optical depths 
Tsph(jc) = T5oo/(x) and rcyi(x) = ts{)i.)J{x), with 



T500 = 1-3530 X 10"= £ 
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and 



C500 - bso^^soo"^ 



Q'yQ'7/(Q'7— ffy) 



(11) 
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Fig. 1. Histogram of predicted values of tsr^,^, the cylindrical 
optical depth times solid angle out io R - 5 Rsoo- 

For the non-isothermal case, we use the average T{x)/f scal- 
ing obtained from the middle panel of figure 3 of Arnaud et al. 
(2010). This scaling is only applied for x < 1, and divides the 
pressure profile to obtain the density, which becomes the inte- 
grand in I{x) and J{x). Since the T{x)/T scaling has only been 
measured for r < R^m {x < 1), a different approach is followed 

for r > RsoQ. In this radius range, we express the electron density 

2/3 

in terms of the pressure and the entropy, K{r) - kT(r)/n^ (r), 
and adopt the relation K{r) cc r'^-^. This defines the scaling of 
density versus radius that enters the outer parts (r > Rsm) of 
the integrals I(x) and J{x). This constitutes our best guess of the 
radial dependence of density in clusters, although in Sect. 5 we 
discuss the motivation and limitations of this approach. Fig. 1 
displays a histogram of the estimated values of the cylindrical 
optical depth integrated out to a radius of R - SR^oo, T5 for 
the non-isothermal case. 

In the analyses described below, we exclude clusters located 
at less than 1 .5 FWHM from point sources detected at more than 
5 cr in any of the single frequency Planck maps. This is done 
in order to remove any spurious signals caused by point sources 
associated with clusters. We also mask clusters lying in regions 
with high Galactic emission, and with estimated masses below 
10'"* Mq. This leads to a basic mask that leaves 1405 clusters on 
the sky (out of the initial 1743 clusters). However, the 2D-ILC 
has its own (and slightly more conservative) mask, which leaves 
only 1321 clusters for analysis on this map. 

In Fig. 2 we show the spatial distribution of the surviving 
clusters. The spatial distribution is quite uniform, except for 
some areas where deeper X-ray observations have allowed for 
more detections. 



2.3. Simulations 

In order to test and disentangle the effects of instrumental noise, 
CMB, Galactic and extragalactic foregrounds on our results, we 
make use of simulations. Specifically, we use 

a) Simulations of clusters with similar characteristics to the 
MCXC sample. We simulate SZ clusters at the actual MCXC 
clusters' locations, using the Amaud et al. (2010) pressure 
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Fig. 2. Locations of MCXC clusters outside the masked region. 
The colour scale indicates inverse noise squared at the posi- 
tions of the MCXC clusters. The large coloured circles indicate 
directions for various dipole determinations: HFI dipole from 
MCXC cluster locations (green); HFI dipole from shifted posi- 
tions (brown); CMB dipole (orange); Kashlinsky et al. (2010) 
dipole (red). The HFI dipole from MCXC cluster locations is 
compatible with CMB and Galactic residuals. 

profile, and the cluster mass and size values obtained from X- 
ray observations. Clusters are assumed to be isothermal. The 
Comptonization Y parameter of each cluster is then com- 
puted using the cluster mass, as outlined in Sect. 2.2. We 
include in the simulations a scatter in the parameters of the 
scaling relations with a normal distribution, which results in 
an averaged scatter in the Y parameter of 10 %. We generate 
a set of 1000 such simulations at the effective frequencies 
in Table 1 in order to assess the effect of that scatter on our 
derived results. 

b) CMB realizations for the current WMAP-1 best-fit model 
(1000 of them). 

c) 1000 noise realizations with noise variance estimated from 
the difference between the first and second halves of the 
Planck rings for a given position of the satellite and divided 
by Planck's appropriate hit map (or exposure map). In do- 
ing so, we disregard noise correlations between pixels. Such 
simulations take into account the non-uniform sky coverage. 
The (non-uniform) noise level in the direction of the MCXC 
clusters is visible in Fig. 2. 

Additionally, we make use of the simulations developed by 
the Planck collaboration (the Planck Sky Model, Delabrouille 
et al. 2012) in order to assess the Galactic contribution to the 
bulk flow measurement on the whole sky. These contain our cur- 
rent best knowledge of the diffuse Galactic component ("PSM 
diffuse"). Since Planck is most sensitive at frequencies above 
100 GHz where the dust emission dominates, we only consider 
this component. We use model number 7 of Finkbeiner et al. 
(1999) which extrapolates the IRAS lOO/zm data to lower fre- 
quency using a modified blackbody frequency dependence and 
a spatially varying dust grain equilibrium temperature based on 
the DIRBE 240/140 jum maps. 

3. Statistical methods 

Here, we outline two different statistical approaches imple- 
mented when searching for the kSZ signal in Planck data. The 



first one (the aperture photometry method) is applied on a clus- 
ter by cluster basis and makes no assumption about the gas dis- 
tribution within a given radius where most of the cluster sig- 
nal must be generated. It nevertheless has to make assumptions 
about the amount of signal that is present in the outskirts. This 
method is quick and easy to implement, in particular when per- 
forming checks for systematic errors. The second approach (the 
unbiased Matched Multi-frequency Filter) makes use of infor- 
mation related to the expected spatial distribution of gas and the 
scale dependence of all sources that can be regarded as noise (in- 
cluding the CMB). The use of these two, different filters is mo- 
tivated by the consistency requirement of having independent 
algorithms that should provide compatible kSZ estimates. The 
unbiased Multi-frequency Matched Filter was implemented in- 
dependently by four different teams, two implementations being 
on square patches centred on clusters, and the other two on the 
whole celestial sphere. The first two implementations target the 
determination of individual cluster peculiar velocities, while the 
latter two are specifically developed to derive local bulk flows, 
since they focus on the dipole pattern of the kSZ effect in clus- 
ters on the celestial sphere. Implementations targeting clusters 
individually allow constraints to be placed on the mean cluster 
velocity (or monopole), the rms or variance of cluster peculiar 
velocities, and also on local bulk flows, i.e., the kSZ dipolar pat- 
tern, as will be shown below. Results from different codes con- 
firm the robustness of our results. 

3. 1 . The aperture photometry method 

The aperture photometry method (hereafter AP) computes the 
average temperature within an input radius R, and subtracts from 
it the average temperature computed in a surrounding ring of 
inner and outer radii R and fR (/ > 1), respectively (see, e.g., 
Hernandez-Monteagudo & Rubino-Martm 2004). In this work 
we use / = ■\/2, so that the outer ring has the same area as 
the inner circle. This is a compromise between having too few 
pixels in thin rings (yielding noisy estimates of the average) and 
being insensitive to local background fluctuations of typical size 
just above R (that are washed out for choices of / which are too 
large). This filter constitutes a simple approach to enhance the 
signal coming from a region of size R against the background. 
In our analyses, the AP procedure was applied in the direction 
of each galaxy cluster, separately in each HFI frequency map. 
When looking at clusters, the filter scale R was taken equal to 
either kO^c^ (k times the angle subtended by the radius at which 
the cluster's density equals 500 times the critical density) or the 
FWHM of the beam, depending on whether the object is resolved 
(A:05O() > FWHM ) or not. Values of k ranging from 0.25 up to 
2 showed that the strongest constraints were obtained for k - 
0.25. Yet smaller values of k do not yield significant differences, 
since for such low k practically all clusters become unresolved. 
We shall describe results with k - 0.25, unless other values are 
explicitly quoted. 

The subtraction of the average temperature in the outer ring 
from the average of the inner circle also removes some fraction 
of the object's flux, which must be accounted for. This results 
in a correction factor of the order of 12 % if clusters are smaller 
than the beam size (FWHM > 6*500). It instead, FWHM < 6*500, 
this correction must make use of some model for the cluster gas 
density profile, as we address next. 

In order to obtain velocity estimates, it is necessary to divide 
the filter outputs (in units of temperature) by the CMB temper- 
ature monopole and a prediction of the clusters' optical depths 
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integrated out to the radial aperture. This then provides an esti- 
mate of the entire cluster's peculiar velocity. The amount of kSZ 
signal that is subtracted by the removal of the outer ring has to 
be accounted for by the same model which provides the opti- 
cal depth versus aperture radius. As explained in more detail in 
Sect. 5, we use the REXCESS observations provided in Arnaud 
et al. (2010) to infer an analytic fit to the average density profile 
in clusters within /?5oo, and use arguments on the behaviour of 
gas entropy at r > /?5oo to extrapolate the density profile at larger 
radii. The adopted model for density provides velocity ampli- 
tudes that are about 28 % higher than those obtained under the 
assumption of isothermal clusters, although we expect clusters 
to be closer to our adopted profile than to an isothermal one. 
Nevertheless, it is our ignorance of the clusters' density profiles 
which drives most of the uncertainties in the velocity constraints. 

When testing for systematic effects, this same filter can eas- 
ily be applied at displaced positions on the sky (that is, positions 
on the sky separated from the real cluster positions by a known 
angle). In the absence of sources and clusters, the average of the 
outputs of this filter at those displaced positions should be com- 
patible with zero, and their rms provides an error estimate for 
the AP filter output at the real cluster's position. 

3.2. The unbiased Multifrequency Matched Filter (uMMF) on 
patches 

The unbiased Multifrequency Matched Filter (Herranz et al. 
2005; Mak et al. 2011, hereafter uMMF) is a Hnear multi- 
frequency filtering technique that is specifically tailored to deal 
with signals that have the same spatial template but different fre- 
quency dependence. A good example of this is the imprint on 
CMB photons caused by the tSZ and kSZ effects. The uMMF 
can be considered as a modification of the Multi-frequency 
Matched Filter (MMF Herranz et al. 2002; Schafer et al. 2006; 
Melin et al. 2006) that optimally enhances one of the two super- 
imposed signals while cancelling out the other. As demonstrated 
in Herranz et al. (2005), it is possible to devise a uMMF that 
detects the tSZ effect and estimates its intensity without the bias 
produced by the kSZ effect, or a different uMMF that extracts 
the kSZ signal and removes the bias caused by the tSZ effect. 
In this paper we are interested in the latter option. In thermody- 
namic units, the uMMF for optimal detection and estimation of 
the kSZ effect is given, in Fourier space, by 



where the constants a, /3 and A are given by 



(12) 
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ay-p^. 
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with y = 









and where P is the cross-power spectrum matrix of the gener- 
alized noise (CMB plus foregrounds plus instrumental noise), 
T - [Ty (k)] is a vector containing the spatial profile of the op- 
tical depth of the cluster (obtained from the universal profile of 
Arnaud et al. (2010) after dividing by the constant temperature 
and convolving by the beam that corresponds to each channel) 
and F - [fvTy (k)] is the vector obtained by multiplying, element 



by element, the profile r by the well-known tSZ frequency de- 
pendence fy. Thus this method observes only isothermal profiles 
for clusters, an assumption which results in a roughly 5 % low 
bias in the radial velocity amplitude, as shown in Sect. 5. The 
power spectrum matrix P is computed from real data in patches 
surrounding the sources. Once we have the filters (Eq. 12), the 
filtered image 



w(x) 



dydi(y)(I>i{x-y) 



(14) 



is optimal for the detection of the kSZ effect and has no trace of 
the tSZ effect. In this equation, c/, represents the unfiltered map 
in the /th frequency channel. The filters are normalized so that 
w{xq), where xo is the location of the centre of the cluster, is 
an unbiased estimator of the kSZ signal due to the cluster An 
estimation of the error of this is given by the square root of the 
variance o-^^,{xo), which can be directly obtained from the filtered 
map or calculated through 



a 
A' 



(15) 



where a and A have the same meaning as in Eq. 13. 

In this work, two different uMMF implementations on square 
patches were used, confirming the robustness of the results. As 
mentioned above, the two implementations assume that the spa- 
tial distribution of the thermal and kinetic signals follows the 
pressure, for which we adopt the universal pressure profile from 
Arnaud et al. (2010). For each cluster, the profile is scaled with 
/?5()0. The two implementations mainly differ in the size of the 
patches used to estimate the background around each cluster and 
the details of the cross-power spectrum estimation on the data. 
In both cases, we apply the resulting uMMF to the patches and 
directly obtain the estimated velocity at the centre of the filtered 
patch. The rms of the filtered patch outside the centre region oc- 
cupied by the cluster gives an estimation of the velocity error 
This leads to a good statistical match between velocities mea- 
sured by the two implementations, but not detailed agreement on 
a cluster by cluster basis. This is expected, since the peculiar ve- 
locity estimate per cluster is dominated by noise, and the actual 
noise component present in each estimate is dependent on the de- 
tails of each specific implementation. The method provides esti- 
mates of the kSZ flux integrated over the cluster profile; these are 
translated into velocity estimates for each cluster after dividing 
by the integrated optical depth. Errors in these estimates of the 
optical depth will lead to errors in the velocity estimates, but, as 
will be discussed below, these should have little impact on esti- 
mates of ensemble quantities like velocity averages, dipoles, and 
rms estimates. More important error offsets (at the 5-25 % level) 
are expected from inaccuracies associated with the gas density 
profile in clusters (see Sect. 5 for details). 

The uMMF method may also be applied to a single map (as 
is the case for the 2D-ILC map), a situation in which the uMMF 
becomes a simple Matched Filter (MF). 

3.3. Constraining kSZ-induced rms in AP/uMMF 
measurements 

Since the signal-to-noise ratio for the kSZ on a typical MCXC 
cluster is very small (see, e.g., (Aghanim et al. 2001) for fore- 
casts on the analysis of bulk flows and the kSZ effect), we at- 
tempt to set constraints on the kSZ signal by performing sta- 
tistical analyses on the entire MCXC cluster sample. We next 
describe our approach to set constraints on the kSZ contribution 
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to the variance of a set of AP/uMMF outputs. This method re- 
lies on the fact that the kSZ contribution to our AP or uMMF 
measurements is uncorrelated with the dominant noise sources 
(CMB residuals, instrumental noise and dust emission). In prac- 
tice this reduces to searching for a kSZ-induced excess variance, 
and this demands a good knowledge of the variance of the vari- 
ance of AP/uMMF measurements, as we next describe. 

In this work, we set constraints on the variance of the clus- 
ter radial peculiar velocities by looking at the variance of our 
filter outputs. For both AP and uMMF filters, the data consist 
of a velocity component (v) plus a number of noise sources 
(CMB anisotropics, instrumental noise. Galactic and extragalac- 
tic emissions not associated with the clusters, etc., here denoted 
by n): 



di - Vi + rii. 



(16) 



where / is the index of each cluster in our sample of N - 1405 
objects (A^ = 1321 under the strict mask). Sufficiently distant 
from such locations, the data d and the noise n coincide. 

We therefore want to measure cr^g^ - Var(v) = 
Tif=i ~ Vm)-^/{N - 1), where v„i is the mean velocity. The vari- 
ance of the data at cluster locations reads; 



= Var(^/) = 



a) 
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Ztli^i " «m)^ EIIi 2(v/ - V,n)(«; - «„) 
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(17) 



where «„, is the mean noise. Assuming that noise terms and clus- 
ter velocities are uncorrelated, for the large number of clusters 
considered here we expect the last term to be subdominant with 
respect to the first two. We therefore interpret the variance of the 
data in Eq. 17 as the sum of variances of the velocity and noise 
terms: 



"-""clusters 



(18) 



We then estimate the noise variance by looking at 100 locations 
near to clusters where noise properties will be similar By do- 
ing so in each of the 100 nearby locations, we can obtain 100 
estimates of the noise variance and hence construct a histogram 
representing its probability distribution. Note that this distribu- 
tion is, in general, not Gaussian. An example from the derived 
noise rms distributions from the AP and uMMF filters are pro- 
vided in the right panels of Figs. 3 and 6, respectively. 

Given the probability distribution of the noise and our mea- 
sured variance at the cluster locations, we can deduce upper lim- 
its for the clusters' velocity variance. Because the variance ve- 
locity term is positive and added in quadrature to the noise, as in 
Eq. 18, we can conclude that at 95 % confidence limit (C.L.) the 
kSZ contribution should be below the following value: 



cr?„(95%) = trL.e.-'^noise(5%) 



(19) 



Here cr^^j^^(5 %) is the noise variance amplitude limiting the 
lowest 5 % of the noise variance distribution. In practice, since 
our histogram is based upon 100 different variance estimates, we 
write: 



<sz(95%) = < 



- cr„. 



,(5th). 



(20) 



In this equation, cr^^^j^^(5th) denotes the fifth lowest AP/uMMF 
output variance estimate picked from the 100 variance estimates 
making the histogram. The quantity crysL^'i' %) constitutes our 



limit of the peculiar radial velocity rms at the 95 % confidence 
level. Such a constraint is, however, obtained from histograms 
built upon only 100 measurements. Using the histograms built 
upon the filter outputs in blank positions we have run Monte 
Carlo simulations and studied the uncertainty on the lower 5 % 
limit on cr^ . if estimated as outlined above. We find that these 

noise 

uncertainties lie typically below the 5 % level when only a sub- 
sample of 100 clusters are used, and below 1 % when using the 
entire cluster sample (around 1400 objects). 



3.4. All-sky bulk flow with the unbiased Multifrequency 
Matched Filter (uMMF) 

In order to evaluate the bulk flow in Planck data, we adopted 
the procedure previously used on simulations for forecasting 
Planck's performance, as detailed in Mak et al. (2011). We 
briefly summarize the approach here, and we refer the reader to 
Mak et al. (201 1) for further information. In this procedure, we 
do not focus on the velocities of individual clusters, but rather 
fit for both amplitude and direction of the bulk velocity for the 
whole ensemble. The first step of the procedure consists of filter- 
ing the observed maps with a whole-sky version of the uMMF 
that adopts the universal pressure profile from Arnaud et al. 
(2010) convolved with the beam profile of a given frequency as 
in the case of uMMF on patches. For this whole-sky version, 
instead of designing the filters individually for each cluster (that 
match its size), we construct one single filter for all clusters, with 
a characteristic scale of 6500 = 8'. This choice is motivated by the 
fact that the average size of the MCXC sample is (dsm) = 7.8'. 
The filtering procedure combines maps at different frequencies 
into a cleaned temperature map that is then used to fit for the 
real spherical harmonic coefficients of the dipole terms (vx, Vy 
and Vz). In doing so, we adopt the effective frequencies listed in 
Table 1. 

We fit the dipole terms of the filtered map using a weighted 
least squares fit that is based on the HEALPix (Gorski et al. 
2005) IDL procedure remove_dipole. We weight the central 
pixels of the clusters that are outside the masked region with 
inverse noise variance weights, i.e., W,- = l/o"Ar,- where cr^,; is 
the ith pixel noise variance calculated from filtered CMB and 
noise realizations. In such realizations, the instrumental noise 
is white and spatially uncorrelated, with a variance estimated 
from the half ring maps and divided by the hit maps appro- 
priate for Planck data in a given pixel. We convert the dipole 
from temperature units {6T) to velocity ones (v) by means of a 
conversion matrix M previously constructed using simulations 
of clusters with the same characteristics as the sample in hand, 
i.e., V = M 6T^. We then evaluate the error on the bulk flow 
dipole coefficients by fitting dipoles to sets of simulations of 
CMB anisotropics, instrumental noise, and the tSZ effect. In or- 
der to do this, we assume that these sources of errors are uncor- 
related, but we consider potential correlations in the errors for 
the dipole coefficients. The magnitude of the dipole velocity fol- 
lows a probability distribution with three degrees of freedom 
that can be computed as follows: 



X 



(21) 



where v is the variable of the distribution, Vm is the mean of 
the velocity as estimated from simulations, and N is the noise 
covariance matrix under consideration. We compute the covari- 
ance matrix by passing 1000 simulations of the noise compo- 
nents (CMB and/or instrumental noise and/or tSZ) through our 
pipeline and performing the dipole fit on them. The scatter in 
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dipole coefficients provides an estimate of the noise correla- 
tions between the dipole directions, i.e., N = ^ynoise^,T,noise^ ^t^^ 
95 % upper limit is then determined to be the velocity at which 

= 7.8, which is the 95 % upper limit forax^ distribution with 
three degrees of freedom. Errors on the bulk flow measurements 
are therefore computed on the basis of simulations, and include 
sources of uncertainties in the mass-observable relation as well 
as in the residual contamination from thermal SZ, CMB and in- 
strumental noise. 

We verified that the most stringent constraints are obtained 
when only the central pixel in the direction of cluster's location 
is considered after filtering the map (and since the data have al- 
ready been matched filtered, applying an aperture would not be 
valid). For the frequency maps used, we present results based on 
the four lowest HFI channels, i.e., 100, 143, 217, and 353 GHz. 
We verified that extending the analysis to the two highest LFI 
channels, 44 and 70 GHz (as in Mak et al. 201 1) gives consistent 
results, but does not significantly improve the constraints. 

Finally, before we end this section we stress the difference 
existing between methods working on patches (such as the ones 
described in previous subsections) and this method, which works 
on the entire celestial sphere. The former methods are insensitive 
to scales larger than the patch size, unlike all-sky methods for 
which filtering is implemented on all angular scales. 



4. Analyses and results 

This section contains the entire set of results in this paper, and is 
divided into three subsections. The first one (Sect. 4. 1) addresses 
the constraints on the kSZ monopole and the rms of the cluster 
peculiar velocities, and is divided into two parts, devoted to the 
results obtained with the AP and uMMF filters. The second sub- 
section (Sect. 4.2) studies the constraints on bulk flows, and is 
divided into four parts. The first three outline the constraints ob- 
tained with the three filters defined in Sect. 3. The fourth part re- 
visits the specific filter implemented by Kashlinsky et al. (2008). 
Finally, the third sub-section (Sect. 4.3) sets constraints on inho- 
mogeneous cosmological models. 

As mentioned above, the MCXC catalogue consists of a sam- 
ple of massive clusters of galaxies hosting large reservoirs of 
hot gas, where the CMB is distorted by means of the tSZ effect 
(Planck Collaboration X 2011; Planck Collaboration XI 2011). 
We target this cluster sample in our attempt to detect or put con- 
straints on peculiar motions in the local Universe. Provided that 
the expected typical correlation length for peculiar velocities is 
of order 20-40/;"' Mpc, we do not expect MCXC clusters to 
show any significant level of coherent motion. Clusters in the 
MCXC catalogue cover a wide redshift range and distances be- 
tween them are far larger than the velocity coherence length. 
However, in the last few years there have been several works 
(Kashlinsky et al. 2008, 2010; Abate & Feldman 2012) claim- 
ing the presence of extremely large-scale bulk flows, and such 
scenarios can be tested with the MCXC cluster sample. 

Note that the linear theory ACDM predictions for the pe- 
culiar velocities of the clusters are supported by the output of 
numerical simulations, although clusters and groups may show 
biases depending on their environment, with higher velocities in 
overdense regions, and non-Maxwellian tails (Sheth & Diaferio 
2001). In any case, the velocity rms expectation, when translated 
into temperature fluctuations via Eq. 1, yields too small a signal 
to be detected on an individual basis (typical velocity estimate 
errors lie at the level of thousands of kms"'). This motivates 



a statistical approach which targets ensemble properties of the 
cluster peculiar velocities. 

We first apply the AP filter to raw HFI frequency maps. Since 
this filter is applied independently on different frequency bands, 
it permits us to track separately the impact of other contribu- 
tions like the tSZ effect or dusty point sources. When impos- 
ing constraints on the cluster peculiar velocities, we also use the 
cleaned 2D-ILC CMB map. Likewise, the use of the uMMF on 
raw HFI frequency maps allows us to test for dust contamination, 
tSZ spectral leakage, or errors in the cluster size determination. 
However, the most restrictive velocity constraints are usually ob- 
tained from the 2D-ILC map. 

4.1. Constraints on kSZ monopole and rms 

In this section we present the constraints that Planck sets on the 
amplitude of the peculiar velocity monopole (average) and rms 
in our cluster sample. 

4.1 .1 . Constraints from the AP filter 

For all MCXC clusters outside the joint HFI mask, an AP es- 
timate is provided for each frequency band. In order to test for 
systematic errors, this filter is applied not only at the cluster po- 
sitions, but also on 100 other positions displaced from the real 
ones in either Galactic or equatorial latitude. For each position, 
the amount of displacement is an integer multiple of three times 
the FWHM of the beam corresponding to the frequency map un- 
der study. For a fixed cluster /, the AP output rms from the dis- 
placed positions provide an estimate of the rms of the AP output 
at the /th cluster's real position {ctap, /). For each displacement, 
we consider only positions outside the effective mask, and com- 
pute both the average (or monopole) of the AP outputs, and their 
rms, as we run over different clusters. The left panel of Fig. 3 
displays the histograms of the AP outputs for the 100 displaced 
positions at 100 GHz (black solid line), 143 GHz (red solid line) 
and 217 GHz (green solid line). The vertical, dot-dashed lines 
correspond to the AP outputs at the real cluster positions (zero 
angular displacement). 

Note that for each displacement some of the real MCXC 
cluster positions may fall in masked pixels. In those cases, the 
AP filter outputs are ignored, that is, for each set of displaced 
positions, the number of usefiil AP estimates equals the num- 
ber of clusters under consideration minus the number of times 
that the "displaced" AP filter centres falls on a masked pixel. We 
hence do not consider AP outputs whenever the filter is centred 
on masked pixels. The left panel of Fig. 3 shows that the AP 
approach is sensitive to the tSZ-induced decrements at 100 and 
143 GHz, since the AP output monopoles at cluster positions fall 
in the negative temperature range, far from the histograms com- 
ing from displaced positions (which are centred near zero). The 
observed monopoles in this panel are less negative (by about 
20 %) than predictions based upon the universal pressure pro- 
file of Arnaud et al. (2010). Given measurement errors, this low 
bias is marginally significant (around 3 cr) and is probably due to 
residual point source emission and/or inaccuracies in the mod- 
elling of the beam impact on our predictions. For 217 GHz, how- 
ever, the AP monopole falls on the positive part of the histogram, 
possibly indicating traces of tSZ-induced emission (since the ef- 
fective frequency of this channel is above the tSZ null). This his- 
togram can be converted into velocity units through dividing by 
each cluster's estimated optical depth (see Sect. 2.2). After av- 
eraging over the full MCXC un-masked cluster sample one can 
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Fig. 3. Colour coding common for all panels; black corresponds to the 100 GHz channel; red to 143 GHz; green to 217 GHz; and blue 
to the 2D-1LC map. Histograms are obtained from AP output at 100 displaced positions, while vertical dot-dashed lines correspond 
to AP outputs obtained on the clusters. The left panel refers to the AP output monopole/average, whereas the middle and right one 
display histograms built upon rms estimates. 



compute the conversion factor from thermodynamic temperature 
fluctuations {6T) to peculiar radial velocity (vy) for this sample, 

v\\^hi.6T. (22) 

This conversion factor however depends on the AP radius 
applied. We obtain values for /t2v of 172kms"^juK"' and 
203kms"' juK"', for the AP radius choices of O.25 05oo and 
0500, respectively. After weighting the AP velocity estimate of 
each cluster by its variance (w, = 1/o"ap ; obtain an es- 
timate for the kSZ average/monopole from HFI 217 GHz data: 
-212 + 80kms"'. If instead we use the 2D-ILC map, the con- 
straint becomes -1 ±73 km s"' . The HFI data in the channel near 
the tSZ null seem to show some residual tSZ contamination (as 
expected from the effective frequency of this channel quoted in 
Table 1), but the 2D-ILC result is consistent with zero. 

The middle and right panels of Fig. 3 display the histograms 
of the rms values obtained from the displaced positions. As for 
the left panel, black and red colours refer to the 100 and 143 GHz 
channels, respectively. In this case, and given the measurement 
uncertainties, predictions from our adopted pressure profile are 
in good agreement with our rms-excess measurements. The AP 
output rms estimates determined at the positions of clusters are 
displayed by the vertical, dot-dashed lines: they fall clearly off 
the histograms obtained from displaced positions, showing an 
excess rms, which is however not seen at 217 GHz (green curve 
in the right panel). In this panel, the AP output rms estimated at 
the cluster positions falls in the middle of the histogram obtained 
from displaced positions. The same occurs for the 2D-ILC map, 
denoted by blue lines. This suggests that the rms excess found in 
the middle panel has come from tSZ, since it does not show up in 
the 217 GHz channel. In ideal conditions, with identical beams 
throughout channels and an absence of noise and foregrounds, 
the histograms in the middle and right panels would be identical; 
the differences among them are reproduced when performing the 
analysis on Planck simulated maps, which account for different 
noise levels and beam sizes. 

Applying the procedure outlined in Sect. 3.3 on the AP rms 
distribution displayed in the right panel of Fig. 3, we set con- 
straints on the kSZ-induced contribution to the total measured 
rms. We find that radial velocity rms constraints for the whole 
cluster sample are, at the 95 % confidence level, 2017 km s"' and 
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1892 km s ', for the raw 217 GHz and 2D-ILC maps, respec- 
tively. These upper limits are about a factor of 8 above theoreti- 
cal predictions, and can be only slightly improved by looking at 
subsets of the cluster sample. We have checked that kSZ errors 
decrease with cluster mass and angular distance, since the more 
massive clusters provide higher kSZ signals and the CMB con- 
tamination is less important on smaller angular scales. However, 
constraints on radial velocities do not improve significantly. 
By using only an un-masked cluster sub-sample containing the 
first 1000 clusters, which have larger values of mass times an- 
gular distance (M500 x Da, <M5oo>subsampie = 2.3 x IO^^Mq, 
(■z)subsanipie — 0.18), the 95% Confidence level constraints 
from HFI 217 GHz and 2D-ILC data become 1806 km s ' and 
1229 km s"', respectively. These are still a factor 5-7 above the- 
oretical predictions for ACDM. 
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Fig. 5. Histogram of recovered radial peculiar velocities as esti- 
mated by the uMMF implementation on HFI frequency maps. 



4.1 .2. Constraints from the uMMF/MF filters 

The uMMF approach can provide accurate kSZ amplitude es- 
timates under the assumption that clusters, as well as being 
isothermal, follow the universal pressure profile of Amaud et al. 
(2010). In Sect. 5 we address the bias that this assumption in- 
troduces when clusters show different density profiles, finding 
a roughly 5 % low bias in the velocity amplitude estimates. On 
an individual basis, the uMMF provides velocity errors that de- 
pend on the mass and the size of each cluster on the sky, and 
lie at the level of a few thousand kms"', with an average value 
of about 4100 km s ' for the un-masked MCXC sample and the 
six HFI channels. When properly accounting for the finite band- 
width of the HFI channel spectral responses, the average pecu- 
liar radial velocity of MCXC clusters is compatible with zero 
(15 + 60kms The 2D-ILC map provides 72 + 60kms 

When binning the cluster sample in redshifts, we again find 
no evidence for any statistically significant average peculiar ve- 
locity (see red and green symbols in Fig. 4). Estimates of the 
kSZ monopole for clusters belonging to different redshift bins 
are shown in Fig. 4. 

Bear in mind that if no colour correction is taken into ac- 
count, and use is made of the nominal HFI frequencies instead 
of the effective ones (see Table 1), then effects associated with 
the finite spectral response in HFI channels become of relevance. 
In Fig. 4 blue squares and black triangles display the average 
radial velocity estimates (within different redshift bins) as in- 
ferred by the uMMF when colour correction is ignored, using the 
three lowest frequency or six frequency channels, respectively. 
In those cases, average velocity estimates lie a few cr below the 
zero level for several redshift bins, pointing to some positive 
residual temperature fluctuations at the cluster positions (which 
is expected since the effective frequency for the third HFI chan- 
nel is about 222 GHz, i.e., above the tSZ null frequency). Colour 
corrections must hence be made when interpreting kSZ measure- 
ments. However, the fact that velocity estimates for three and six 
HFI bands are compatible suggests that dust contamination is 
not important. We also checked that no significant differences 



are found when introducing changes of + 1 cr in the adopted HFI 
Gaussian fit FWHM values. 

In Fig. 5 we display the histogram of radial peculiar veloci- 
ties as estimated, for MCXC clusters, by the uMMF implemen- 
tation. The distribution is almost symmetric around zero and has 
a standard deviation of 4100 km s"' . It shows non-Gaussian tails, 
presumably due to the presence of un-resolved point sources and 
other non-Gaussian signals in some galaxy clusters. 

We next conduct a spatial analysis of the Matched Filter 
outputs, just as done for the AP approach above, and again 
analysing both the set of HFI frequency maps and the 2D-ILC 
map. In this analysis we first apply the filter on 100 blank, ran- 
dom positions on the sky where no MCXC clusters are found, 
and compare the resulting rms of the filter velocity estimates 
with the rms of the estimates obtained for the real cluster posi- 
tions. Contrary to the displaced positions for the AP filter, in this 
case the random positions on which the uMMF/MF filters are 
evaluated are not displaced in either equatorial or ecliptic lati- 
tude with respect to the real cluster, but simply randomly placed 
within a 10° x 10° patch centred on the real object. If the kSZ 
signal generated for those "sources" leaves a measurable effect, 
then the rms at the real positions of the clusters must be larger 
than the rms obtained at blank positions. 

In each patch centred on each MCXC cluster, we record 
the uMMF/MF outputs for 100 random positions not coinci- 
dent with the centre. This provides peculiar velocity estimates 
for positions where we expect the kSZ effect to be zero while, 
at the same time, having similar levels of instrumental noise and 
foreground contamination as the positions corresponding to real 
MCXC clusters. For each set of 100 positions, it is possible to 
compute the rms and compare it to the expected value predicted 
by the uMMF/MF method, as provided by Eq. 15. This is shown 
in the left panel of Fig. 6 for HFI frequency maps: the solid line 
displays a one to one relation, and it is roughly followed by 
the recovered rms from the random estimates within the patch 
(vertical axis) versus the predicted velocity errors (horizontal 
axis). Here we neglect all uMMF/MF outputs within the patch 
that fall in masked pixels. After fixing one of the 101 positions 
within each patch, we compute an average velocity by consid- 
ering velocity estimates in all patches at that particular position. 
The histogram of these average velocities computed in the 100 
displaced positions is shown in the middle panel. The vertical, 
red dashed line corresponds to the average velocity as computed 
from the velocity estimates at the patch centres, that is, at the 
real MCXC cluster positions. As expected, the middle panel of 
Fig. 6 shows that the entire ensemble of MCXC clusters exhibits 
average peculiar velocities that are compatible with zero. The 
right panel displays the histogram of the velocity rms estimates, 
computed exactly in the same way as for velocity averages. For 
each of the 101 positions within a patch, we calculate a veloc- 
ity rms after considering the filter outputs for the whole set of 
patches at that position. We end up with 100 velocity rms esti- 
mates in displaced positions, and the resulting histogram in the 
right panel is compared to the velocity rms estimated at the clus- 
ter positions, again shown as a vertical, red dashed line. From 
this distribution, and after following Eq. 20, the uMMF/MF can 
set the following upper limits at 95 % C.L. on the cluster radial 
peculiar velocity rms: 1514 km s"' for HFI frequency maps; and 
987 km s ' for ILC data. As for the AP filter, when restricting 
ourselves to the cluster sub-sample which maximizes the prod- 
uct M500 xDa, we obtain improvements on these constraints: for 
the 1000 clusters maximizing M500 x Da, upper limits become 
798 and 754 km s ' for the raw HFI and 2D-ILC data, respec- 
tively. If we instead choose the top 100 clusters in the sorted 
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Fig. 6. Left panel: Estimated rms of uMMF outputs at 100 random positions within the patch surrounding each MCXC cluster 
versus the corresponding peculiar velocity error as predicted by the uMMF. Middle panel: Histogram of mean kSZ radial velocity 
estimated using all un-masked MCXC clusters, for the 101 positions considered in each cluster patch. The vertical, red dashed line 
corresponds to the average kSZ estimate when averaging throughout the patch centres, i.e., at the real cluster positions. Right panel: 
uMMF velocity rms histogram obtained for the MCXC cluster set, at different positions within the patches, just as for the middle 
panel. Again, the vertical red dashed line corresponds to the kSZ velocity rms estimated at real cluster positions (on patch centres). 



list of the previous sub-sample, the constraints change, but only 
slightly: 794 and 614km s ' for raw HFl and 2D-1LC maps, re- 
spectively. These limits have systematic uncertainties at the level 
of a few percent, but are nevertheless a factor of about 3 above 
the level of ACDM predictions. 



4.2. Constraints on bulk flows and the kSZdipole 

In this section we describe the constraints that Planck sets on the 
existence of bulk flows at different scales and the measurement 
of the kSZ dipole in our cluster sample. 



4.2.1. Constraints from individual cluster velocities 

Extensive efforts have been made in recent years to try to set 
constraints on the local bulk flow (Hudson et al. 2004; Watkins 
et al. 2009; Feldman et al. 2010; Nusser & Davis 2011; Ma & 
Scott 2012), without reaching full agreement so far The kSZ es- 
timates from Planck provide a different approach to the question 
of the local bulk flow: if clusters embedded in structure around 
the Local Group are comoving with it towards a nearby over- 
density, then the kSZ measured for those sources should show a 
dipolar pattern. By looking at clusters within spheres of differ- 
ent radii from us it is possible to set constraints on the local bulk 
flow within different volumes. This provides a direct test on the 
studies of Kashlinsky et al. (2008, 2010, 2012), which claim that 
clusters extending at least up to 800 /i"' Mpc are part of a bulk 
flow of amplitude about 1000 km s"' . 

Given any MCXC cluster sub-sample, we compute the am- 
plitude of the kSZ dipole along a given direction n dip by mini- 
mizing = Yjj{v\\j - a(Mdip ■ n j))' jcrl^. Here, \\j denotes the 
AP/uMMF/MF radial velocity estimate of the y'th cluster (which 
is located in the direction h j), and a-^. is its associated error 
After assuming uncorrelated errors (from cluster to cluster), this 
minimization yields both an estimate of a and a formal associ- 
ated error: 
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Fig. 7. Upper limits at the 95 % confidence level for the dipole 
amplitude from MCXC clusters contained in local spheres of 
varying radii. Blue and black arrows denote hmits for AP and 
uMMF methods, respectively. The upper hmits are indicated by 
the tails of the arrows. 



cr,-v - 



2/«dip ■ h jYla-i,^ 



(23) 



We compute the kSZ dipole for cluster sub-samples contained 
within increasing radii from the Local Group. For each cluster 
sub-sample and associated kSZ estimates, we fit a dipole along 
every direction in the sky, i.e., we sweep in Mdip, and retain the 
direction which yields the highest uncertainty (i.e., highest cr,-, 
value) in order to find the corresponding upper limit. Figure 7 
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Fig. 8. Mollweide projection in Galactic coordinates of the upper 
limit (at 95 % C.L.) of the kSZ dipole amplitude from applying 
the uMMF approach to HFI frequency maps using the whole 
MCXC cluster sample. In no direction is the dipole detected at 
more than 2 cr. 



displays the corresponding upper Umits (at the 95 % confidence 
level, calculated as 2cr, assuming a Gaussian distribution) on the 
dipole values for different radii and both AP and uMMF meth- 
ods on HFI frequency maps: for both methods, the 95 % C.L. 
upper limit for radii of 90 /i ' Mpc amounts to ~2000kms"', 
but it decreases rapidly as the volume increases. For spheres 
of radius around 350 /z"' Mpc, the uMMF limits fall to about 
390 km s"' , and the corresponding AP limit is just slightly higher 
(520km s"'). In the largest volume probed by the MCXC clus- 
ters, the 95 % C.L. upper limits become 329 and 254 km s"' for 
the AP and uMMF filters, respectively. Despite being very differ- 
ent in their definition, the two methods give rise to a very similar 
pattern in the bulk flow constraints inside different volumes, and 
in all cases the measured dipoles are compatible with zero. 

In Fig. 8 we display the 95 % upper limit on the kSZ ampli- 
tude from the uMMF filter using the whole MCXC cluster set 
(for which {z) - 0.18) and HFI frequency maps. In no direction 
does the measured dipole exceed 2cr, and the direction with the 
highest ajcTe, value is close to the Galactic plane. This is to be 
expected if the errors in the Galactic x- and y-components of the 
dipole are larger than the z-component, due to the lack of clusters 
at low Galactic latitudes. When restricting ourselves to clusters 
below z - 0.25, the dipole amplitude along the CMB dipole 
direction (/,/?) = (264°, 48°) (Hinshaw et al. 2009) amounts to 
80+ 150 km s"', and limits to 50 ± 160 km s"' along the direction 
of apparent motion of the Local Group with respect to the CMB, 
iUb) = (276°,30°)(Kogutetal. 1993). This result is in clear con- 
tradiction with Kashlinsky et al. (2010), who find a bulk flow of 
about 1000 km s"' amplitude within radii of 300-800 /r' Mpc at 
~ 3 cr C.L. Since our error bars are a factor of about 2 smaller, 
this suggests that we test the outcome of the filter used by those 
authors on our data (see Sect. 4.2.3 below). 

4.2.2. Constraints from all-sky method 

We can also compute the bulk flow according to the procedure 
outlined in Sect. 3.4. We filter the observed HFI maps and fit 
monopole and dipole velocity coefficients to the filtered data, as 
well as to simulations of the data (PSM diffuse component, tSZ 
and CMB plus instrumental noise simulations; see Sect. 2.3). 
Simulations provide one way of estimating uncertainties which 



allow one to propagate errors on cluster-derived measurements, 
mainly induced by dispersion in the scaling relations, throughout 
the whole pipeline. An alternative procedure to derive uncertain- 
ties consists of randomizing the positions of the clusters on the 
sky and computing the monopole and dipole from these random 
directions, adopting the same procedure used on the real data. 
We use the latter to show the typical variations of the diffuse 
component's contribution for small displacements around clus- 
ter locations. Specifically, we consider directions displaced by 
30' to 1° from the cluster nominal locations, while also avoiding 
mask boundaries (these are the "shifted positions" in Table 2). 

Values for the resulting velocity dipole coefficients are pre- 
sented in Table 2. The main result is that Planck data give dipole 
coefficient amplitudes consistent with those expected from the 
ACDM scenario, once one has taken into account the contamina- 
tion from Galactic foregrounds and other signals. The apparent 
bulk flow measured is 614 km s"'. However, with this particular 
configuration for cluster positions, the diffuse Galactic compo- 
nent provides a non-negligible contribution to the dipole signal, 
529 km s ', as measured in the PSM simulations. The errors on 
the diffuse component, as estimated by randomizing the clus- 
ter directions on the PSM diffuse component simulations, are 
smaller than those induced by the thermal SZ and CMB plus 
instrumental noise simulations (see Table 2). 

Simulations of the tSZ component, which account for un- 
certainties in the SZ signal for clusters with a given temperature, 
induce a 1 cr uncertainty on the bulk velocity of 40 km s" ' , and an 
overall bias in the velocity estimation of the order of 400 km s ' . 

Uncertainties from CMB confusion and instrumental noise 
( 140-290 km s ' in the different directions) are dominant over 
tSZ ones. The fraction of the observed bulk flow not accounted 
for by Galactic foregrounds (by subtracting the dipole as a vec- 
tor, this amounts to 350 km s"') is within 95 % of the error on 
bulk flows induced by the tSZ, CMB and instrumental noise 
(893 km s ') and below the 95 % level of CMB plus instrumental 
noise alone (543 km s"'). 

By restricting the cluster sample to the objects within a spec- 
ified distance from us, it is possible to constrain the bulk flow 
within spheres of a given comoving radius. This is what is dis- 
played in Fig. 9, where the Galactic component has been sub- 
tracted. We notice that, for all distances, the measured bulk 
flow is below the 95 % confidence level as measured from 
maps including only CMB, instrumental noise, and tSZ clusters. 
The upper limits reach an approximately constant value above 
scales around 500/;"' Mpc, as a small fraction of the clusters 
in this sample are at larger distances. The 95 % upper limits at 
2400/2 ' Mpc are 893 km s"' when all sources of noise are con- 
sidered, reducing to 543 km s ' when CMB plus instrumental 
noise are taken into account. 

The results reported in Fig. 9 refer to the nominal mask, 
while in Table 2 we also quote results for the more restrictive 
mask. The two sets of results are very similar, however 

In this analysis, we also fit for the direction of the measured 
bulk flow. Even although the detection is not significant, it might 
still be instructive to compare the best fit direction to other po- 
tentially relevant directions. Results for various cluster configu- 
rations and Planck data are displayed in Fig. 2, together with the 
CMB dipole and the claimed dipole direction of Kashlinsky et al. 
(2008). We notice that the direction we determine from Planck 
data and MCXC clusters is quite different from both the CMB 
dipole and the result of Kashlinsky et al. (2008). It aligns better 
with the direction of the collection of clusters in the map, which 
happen to be in a low instrumental noise area of the sky, as one 
would expect from a noise-induced measurement. Indeed, sim- 
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Fig. 9. Bulk flow amplitude measured in Planck data with the all- 
sky method, after subtraction (vectorially) of the Galactic contri- 
bution (black crosses), compared with 95 % upper limits derived 
from simulations containing CMB and instrumental noise only 
(blue arrows) or also including tSZ signal (black arrows). The 
fact that the crosses are below the arrows at all scales shows that 
there is no significant bulk flow detection. 



ulations show that the directions of bulk flows of the magnitude 
seen in the data cannot be recovered with great precision. Errors 
are of the order of tens of degrees, depending on the bulk flow 
direction (Mak et al. 201 1). 

Finally, we notice that the upper limits to the bulk flow that 
we find with this method are above those found in the previ- 
ous subsection. This is not surprising, as we are fitting here for 
both the velocity direction and amplitude, and we compute er- 
rors in a different way. The upper limits obtained with this ap- 
proach should be considered as more conservative. Nevertheless 
they are about a factor of five better than what was found using 
WMAP data. 

4.2.3. Revisiting tine Kaslilinsky et al. (2010) filter 

The idea of constraining the local bulk flow of matter by look- 
ing at the dipolar pattern of the kSZ in the galaxy cluster pop- 
ulation was first discussed by Haehnelt & Tegmark (1996) and 
further developed by Kashlinsky & Atrio-Barandela (2000). The 
method was applied by Kashlinsky et al. (2008, 2009) to WMAP 
data, analyses that have been followed by more recent studies 
(Kashlinsky et al. 2010, 201 1). In this section, we perform a di- 
rect application of their filter to both WMAP and Planck data, 
and interpret it at the light of the results already outlined in this 
work. 

We first implement the filter of Kashlinsky et al. (2010) on 
the MCXC cluster sample and the WMAP-7 data. After using the 
extended temperature KQ75 mask, we obtain filtered maps from 
the cleaned Q, V and W band WMAP data. Since the filtered 
maps for the four W-band Differencing Assemblies (DAs) used 
by those authors are publicly available^^, a direct comparison of 

^ The data were downloaded from the URL site 
http : //www . kashlinsky . info . 



the filtered maps can be performed: for instance, for the filtered 
maps corresponding to the fourth W-band DA, the temperature 
rms outside the joint mask in our filtered map is 74 //K, very 
close to the 77 fiK obtained from the map used by Kashlinsky 
et al. (2010). The rms of the difference map amounts to 35 juK, 
and a visual inspection shows the similarity between both maps. 
Each cluster is assigned a radius of 25', and the remove_dipole 
routine from HEALPix is used when computing the monopole 
and dipole in the subset of pixels surrounding the clusters. The 
monopole and dipole components obtained for the WMAP W 
band are displayed by the black, vertical dot-dashed lines in 
Fig. 10. These are in very good agreement with the results ob- 
tained by Kashlinsky et al. (2010). 

We next distribute the same number of clusters surviving the 
mask randomly on the unmasked sky 1000 times, assign them a 
circle of radius 25' and repeat the monopole and dipole compu- 
tation. For each of the 1000 cluster configurations, we separately 
compute the monopole and dipole for each of the DAs. This per- 
mits us to obtain the rms for each component and DA, in such 
a way that a combined estimate of the monopole and dipole can 
be extracted from all DAs by inverse-variance weighting the es- 
timate for each DA. This is caiiied out for the real cluster con- 
figuration on the sky and for the 1000 mock (random) configura- 
tions. From the latter, we obtain the histograms shown in Fig. 10. 
The average quantities out of the 1000 simulations are displayed 
by the solid, vertical lines. Black lines refer to WMAP data, and 
our results show that the y-component of the dipole is peculiar, 
in the sense that it falls far in the negative tail of the distribution. 

When repeating these analyses with the 2D-1LC map, we ob- 
tain the results displayed by the red lines in the same figure. In 
this case, the dipole components from the real data fall further 
outside the distribution provided by the histograms, as none of 
the 1000 mock cluster configurations provides a dipole of larger 
amplitude than the one measured from the real MCXC sample. 
These results suggest that the dipole measured at the MCXC 
cluster positions is indeed peculiar if compared to dipole esti- 
mates from randomized cluster positions. 

Nevertheless, there is one aspect to be studied more closely, 
namely the angular distribution of clusters on the sky. In what 
follows, the filtered map built upon the 2D-1LC data will be used. 
So far our Monte Carlo simulations assumed that clusters were 
placed randomly on the sky, i.e., the clustering of our sources 
has been neglected. We next perform tests in which the angu- 
lar configuration of our MCXC cluster sample is preserved. The 
first test consists of repeating the filtering and subsequent dipole 
computation on 1000 CMB mock skies following the WMAP- 
7 best-fit model. These mock CMB maps contain no kSZ and 
hence should give rise to no significant kSZ dipole. Out of this 
ensemble of mock skies, we compute the dipole using the posi- 
tions of MCXC clusters (as described above) and obtain a his- 
togram from the recovered dipole amplitude. This permits us to 
judge how peculiar our measurement is with respect to the sim- 
ulation outputs. In a second test, we rotate the clusters' angular 
positions around the Galactic z-axis on the real filtered map ob- 
tained from the 2D-ILC data. We conduct 360 rotations of one 
degree step size, in Galactic longitude, while preserving Galactic 
latitude, and the relative angular configuration of MCXC clus- 
ters on the sky. Since the mask mostly discards pixels at low 
Galactic latitude (close to the Galactic plane), most clusters that 
are originally outside the mask remain outside the mask after ro- 
tating. For each rotation a new value of the dipole is recorded, 
and information on dipole statistics is then built up using outputs 
obtained from the real map with the real rotated cluster config- 
uration on the sky. This rotation test, unlike the one based upon 
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Fig. 10. Monopole and dipole component estimates after applying the spatial filter of Kashlinsky et al. (2008) on WMAP-1 W-band 
data (black lines) and on Planck 2D-ILC maps (red lines). Estimates from real MCXC clusters are displayed by vertical, dot-dashed 
lines. Histograms are obtained after repeating the analysis on 1000 random cluster configurations, with averages indicated by the 
vertical, solid lines. The y-component appears discrepant here, but compare Fig. 1 1 . 



CMB mock skies, accounts for the impact of noise, foregrounds 
and other systematic signals that may be present in the filtered 
map. 

The results are shown in Fig. 11. The black histogram re- 
flects the statistics of the recovered dipole amplitudes after 
drawing 1000 random cluster configurations on the real filtered 
map, just as done for Fig. 10. The blue, triple dot-dashed his- 
togram corresponds to the dipole outputs obtained after rotating 
in Galactic longitude one single random cluster configuration ap- 
plied to the filtered map obtained from 2D-1LC data. Clearly, this 
rotation gives rise to a histogram that is very close to the one 
obtained from the 1000 random cluster configurations. On the 
other hand, the green, dot-dashed histogram reflects the statis- 
tics of the recovered dipole amplitudes obtained from the 1000 
Monte Carlo CMB simulations. Again, this histogram is fairly 
close to the one obtained after rotating the real cluster sample 
in Galactic longitude on the real filtered map (red dashed his- 
togram). The recovered dipole amplitude from the real cluster 
positions on the real filtered map is displayed by a vertical black 
line. 

While the measured dipole falls in the far positive tail for 
the simulations using random cluster configurations (black and 
blue histograms), it is however quite unremarkable when com- 
pared to the simulations accounting for the real configuration 
of clusters on the sky (red and green histograms): about 1 1 % 
of cases in both CMB simulations and rotations yield dipoles 
larger in amplitude than the one measured on the real data. The 
green histogram shows that the apparent dipole can be explained 
by chance alignments of random, uncorrected CMB skies. The 
impact of instrumental noise and other component only shifts 
the histogram slightly (as a comparison of the red and green 
histograms suggests). These results show that the dipole mea- 
sured for the real MCXC cluster positions is not peculiar when 
compared to other dipole computations, either on mock CMB 
skies or on the filtered 2D-1LC map for a set of positions in 
which the angular clustering of the MCXC sample is preserved. 
When repeating this analysis on a sub-sample of MCXC clus- 
ters containing the 200 most massive objects, wider histograms 
from both rotations and CMB mock skies are obtained. The area 
under the histograms above the apparent dipole obtained from 
real data at zero-lag rotation amounts to about 56 % of the to- 
tal, see Fig. 12. Unlike for the entire MCXC cluster sample, the 
histograms obtained after randomly distributing this sub- sample 
of massive clusters on the un-masked filtered map are very sim- 



ilar to those obtained after running CMB mocks or rotating the 
clusters in Galactic longitude. This is in better agreement with 
Atrio-Barandela et al. (2010), who found no significant differ- 
ence between the histograms obtained from CMB mocks and 
from randomly distributing clusters on the filtered map. This is 
likely due to the absence of any significant intrinsic dipole in the 
angular distribution of this (smaller) cluster sub-sample. 

Finally, we perform a direct comparison of our results with 
Kashlinsky et al. (201 1). For this purpose, we use the WMAP fil- 
tered maps and the sky mask for galaxy clusters at z < 0.25 used 
by Kashlinsky et al. (2008); these data are presented as supple- 
mentary materials for Kashlinsky et al. (2011). For the filtered 
maps corresponding to the four DAs of the WMAP W band data, 
we apply the rotation test in Galactic longitude. We find that, al- 
though the y-component of the dipole at no rotation is marginally 
peculiar (at the roughly 1-3 % level), the amplitude of the dipole 
is not, since around 14 % of the rotations yield higher amplitude 
dipoles. This is in good agreement with the Planck results out- 
lined above. Hence, according to our estimations of the dipole 
uncertainty, we conclude that the roughly 3 fiK dipole measured 
by Kashlinsky et al. (2008) should not be assigned to the clus- 
ters' peculiar motion, but rather to residuals (mostly of CMB 
origin) in the filtered map. 

4.3. Constraints on inhomogeneous cosmological models 

The sensitivity of the kSZ effect to peculiar velocities and 
bulk flows makes it an excellent probe of nonstandard inhomo- 
geneous cosmological models (Goodman 1995). In particular, 
models in which we are located near the centre of a spherically 
symmetric Hubble-scale void have been examined extensively 
in recent years as alternatives to standard accelerating mod- 
els on Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
grounds (see, e.g., Clarkson 2012, and references therein). Such 
void models can easily reproduce the Type la supernova (SN) lu- 
minosity distance-redshift data without dark energy or modified 
gravity. However, these models generically predict very strong 
outwards-directed bulk flows, due to the greater expansion rate 
within the void, and so are expected to produce a large kSZ 
monopole signal (superimposed, of course, on the usual kSZ 
signal from structure). Constraints on such models using kSZ 
upper limits from nine clusters (Holzapfel et al. 1997; Benson 
et al. 2003; Kitayama et al. 2004) show that the largest voids 
are at odds with these early data, assuming purely adiabatic ini- 
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Fig. 11. Histograms of the recovered cluster dipole amplitude: 
(a) from the 1000 Monte Carlo random cluster configurations on 
the sky (black, solid line); (b) from rotating one random cluster 
configuration in Galactic latitude on the real filtered map (blue, 
triple dot-dashed line); (c) from rotating the real MCXC cluster 
configuration around the Galactic z axis on the real filtered map 
(red, dashed line); and (d) from applying the filter on the position 
of our MCXC cluster sample in 1000 Monte Carlo CMB simu- 
lations following the WMAP-1 best-fit model (green, dot-dashed 
lines). The dipole amplitude recovered at the real MCXC cluster 
positions on the real filtered map is shown by the vertical, solid 
line. It is not significantly detected, provided one is careful to 
simulate the most important noise contributions. 
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Fig. 12. Same as Fig. 11, but restricting the analysis to the 200 
most massive clusters outside the mask. The colour and line- 
style coding is identical to that figure. In this case, all histograms 
contain the measurement on the real positions of clusters (dis- 
played by the vertical solid line). 



tial conditions (Garcfa-Bellido & Haugb0lle 2008; Yoo et al. 
2010). Tight constraints have also been imposed using upper 
limits on the kSZ power from small-scale CMB experiments 
(Zhang & Stebbins 2011; Zibin & Moss 2011). However, the 
results based on small-scale kSZ power are uncertain, due to our 
inability to properly perform perturbation theory in void models 
and our lack of knowledge about the small-scale matter power 
spectrum and baryonic physics (Zibin & Moss 2011). The very 
tight Planck constraints on the kSZ monopole presented above 
are therefore expected to provide extremely stringent limits on 
any such large-scale features, in a manner that is free of the un- 
certainties due to small-scale structure. 

We first briefly describe our void models and calculation 
methods; full details can be found in Moss et al. (2011). 
Growing-mode void models are characterized by a single ra- 
dial function, e.g. the matter density profile. Models with sig- 
nificant decaying modes are ruled out by their extremely large 
kSZ and CMB spectral distortions (Zibin 201 1; Bufl et al. 2012) 
and so will be ignored here. In this study we consider a family 
of smooth void profiles (taken from Moss et al. 201 1) parameter- 
ized by a width, L, and a depth, 6q < 0. Explicitly, we superpose 
(at early times) on a spatially flat background the total matter 
density contrast profile 




r<L, 
r> L, 



(24) 



for comoving radial coordinate r centred on us. In order to ex- 
press our constraints in terms of more directly observable quan- 
tities, in place of L we use the corresponding redshift, zl, at 
which we observe an object at r = L. In place of the depth, 
6o, we use the local matter density parameter at the origin today, 
= 8;tGpm,o/(3//q), where pM,o is the current total matter 
energy density at the centre. Thus Q}^ generalizes the familiar 
density parameter of an FLRW cosmology, and the deepest voids 
have the smallest values of fl}^'^. 

These models require a relativistic treatment at late times, 
since they are not well described by small perturbations from an 
FLRW background. We use the Lemaitre-Tolman-Bondi (LTB) 
exact solution to Einstein's equations to calculate the radial ve- 
locity, V||_ltb(z), between a comoving scatterer (i.e., a cluster) and 
the local CMB frame, using the method of Moss et al. (201 1). 

We compare our LTB void models with the AP Planck radial 
velocity estimates for 1405 clusters by calculating the likelihood 
while varying the width and depth of the void profile. The like- 
lihood, X., is calculated using the full error distributions from 
the 100 displaced positions for each cluster Importantly, for the 
deepest and widest void models, the typical velocities V||,ltb are 
outside of the range of velocities for the displaced positions. This 
means that we cannot calculate the actual likelihoods for these 
models (because they are so small). Instead, we conservatively 
assign such models likelihoods based on the outermost sampled 
regions of the error distributions. This will almost certainly re- 
sult in a large overestimate of the likelihood for these models, 
and hence we will underestimate the confidence at which they 
are ruled out. 

In Fig. 13 we plot contours for the quantity logio(X/-Chom) 
in the width-depth plane. Here -Chom is the likelihood for the 
exactly homogeneous model, i.e., the model for which <5o = 
(which implies Q}^ - 1 and vy ltb(z) = 0). Also shown on the 
plot are the confidence levels for the same void models, but us- 
ing the Union2 compilation of Type la SNe, taken from Zibin & 
Moss (201 1). The SN data demand deep (i.e., low D}^) but wide 
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Fig. 13. Solid contours indicate \og\Q(JLI JLhom) for the AP 
Planck frequency maps as a function of central matter density 
parameter, and the width of the void in redshift, Zh- The 
deepest voids have the smallest values of Qjj^'^. Dotted contours 
are the 1, 2, and 3cr confidence levels from Zibin & Moss (2011), 
using the Union2 SN data. Void models which fit the SN data are 
ruled out at very high confidence by the kSZ data. 




Fig. 14. The solid, red line displays the ratio of kSZ fluxes com- 
puted in circles centred on clusters with radii V2 6500 and 6500 
for different scalings of the entropy with radius, K{r) cc . The 
green, dashed line observes the isothermal case, in which den- 
sity exactly traces the universal pressure profile adopted in this 
work. 



void profiles, while the PlanckkSZ data rules out all but the very 
shallowest (i.e., Q}^ ^ 1) or narrowest (i.e., small zl) profiles. 
Adiabatic void models are thus ruled out at extremely high con- 
fidence. It is easy to understand the strength of this result: voids 
fitting the SNe have V||,ltb(z) ~ 10"* km s"' at z ^ 0.5 (Moss et al. 
201 1), which places them a few standard deviations into the tails 
of roughly 1000 cluster measurements. 

5. Robustness of the results 

Returning to our main goal of determining the kSZ peculiar ve- 
locity constraints, we now address the sensitivity of our results 
to uncertainties in the density profile adopted for clusters and 
the errors in estimates of the optical depth. We also quantify the 
impact of non-CMB noise sources on our error bars. 

5.1. Impact of changes in the density profile 

As mentioned in Sect. 2.2, the adopted density profile at radii 
below /?5()o is a fit from REXCESS observations, but at larger 
radii (x > 1) the density is expressed in terms of the gas pressure 
(which follows the universal profile of Arnaud et al. 2010) and 
the entropy K(r) = kT(r)/nl^^(r). As well as being physically 
motivated, the reason for this is the existence of some constraints 
on the scaling of the entropy for r > R500. As long as the shock 
front is beyond 5 x Rsqq (i.e., beyond our upper radial integra- 
tion limit), the entropy should increase with radius with a power 
law of the type K{r) cc , with y = 1 . 1 for the adiabatic case, 
as predicted by Voit et al. (2005). Observations at r < /?5oo are 
reasonably close to this prediction, but at larger radii, the scaling 
should become shallower (as suggested by the observations of 
Walker et al. 2012). Our computations for the AP filter adopted 
the fixed value y - 0.5 for r > but we explored the impact 
of different scalings when estimating the ratio of fluxes inside 
the inner and outer circles of the AP filter. The red solid line in 



Fig. 14 displays the ratio of the kSZ flux computed in circles 
of radii V2 6*500 and 6500 for different scalings of the entropy, 
K(r) oc r''. This ratio depends on the profile, but not on the clus- 
ter mass or redshift. The horizontal, dashed green line displays 
the isothermal case at all radii, including for x < I. When trans- 
lating these flux ratios into velocity constraints, we find that the 
y = 0.5 profile introduces a boost in the pecuUar velocity am- 
plitude estimate of about 28 % with respect to the case where 
clusters are assumed to be isothermal. For the extreme case of 
y = this error amounts to 36 %, while for y - 1 . 1 it goes down 
to 20%. With respect to our reference model of y = 0.5, we 
expect that variations of y in the range [-1.1,0] will introduce 
changes in the velocity constraints at the level of +6%. If the 
pressure profile is changed to use the parameters that best fit the 
external profiles measured out to 3 x R^oo (Planck Collaboration 
V 2012), then the change in the velocity constraints are small 
enough to be considered negligible. 

Regarding the uMMF, only the isothermal profile was im- 
plemented in the filter. In order to assess the impact of this ap- 
proximation, we conduct a numerical experiment, consisting of 
assigning a 1000 km s"' amplitude dipole to the cluster set of 
our simulations, and extracting this dipole after injecting both 
the isothermal and non-isothermal profiles to the clusters, but 
considering only the isothermal profile in the uMMF defini- 
tion. We recover amplitude values of 1016 ± 106kms ' and 
966 + 106kms ' for the isothermal and non-isothermal pro- 
files, respectively. The errors do not change, which is expected, 
since the profile coded in the uMMF is isothermal in both cases. 
The difference in the recovered values, however, suggest a po- 
tential bias in the amplitude estimate of about 50kms"', which 
amounts to around 5 % of the signal amplitude. This test in- 
dicates that the uMMF is less sensitive to the assumed den- 
sity profile than the AP, probably because the former assigns 
more weight to the central regions of the cluster where the noise 
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Table 3. Impact of uncertainties in estimates of the optical depth 
of clusters on the constraints imposed on the kSZ monopole, 
rms, and dipole from Planck data. Different t error amplitudes 
(given by cr^ ) and two cluster (sub-)samples are considered here. 
The percent levels correspond to the fractional changes found in 
the kSZ monopole error bar (third column), the 95 % confidence 
limit for the kSZ-induced rms-excess (fourth column), and the 
kSZ dipole error bar (fifth column). 
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0.4 
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(mostly induced by the CMB) is lower and where the decreasing 
radial pattern of temperature is less dramatic. 

Hence, we conclude that uncertainties in the radial gas dis- 
tribution in galaxy clusters should introduce errors in the Umits 
imposed on peculiar velocities at the 5-30 % level. 

5.2. Other sources of error 

After repeating the analysis of the AP and MF filters on a pure 
CMB simulated map, we found that errors decreased to about 
70 % of their amplitude on HFI raw maps. This shows that the 
main limiting factor when estimating kSZ velocities is the in- 
trinsic CMB component. The presence of point sources, instru- 
mental noise and other foregrounds should be included in the 
remaining roughly 30%, and after properly accounting for the 
spectral response of HFI detectors, there seems not to be any 
significant tSZ leakage biasing the peculiar velocity estimates. 

The AP and uMMF/MF velocity estimates, however, rely 
on an accurate knowledge of the cluster optical depths within 
a given radius. For each cluster, the integrated optical depth is 
estimated from the adopted radial density profile, which itself 
relies on the Fsoo-A^soo and T-Msm relations. Apart from the 
uncertainties in the shape of the profile (addressed above), the in- 
trinsic scatter in these scaling relations could have an impact on 
the peculiar velocity estimates. In order to test this we conduct a 
Monte Carlo analysis consisting of introducing un-correlated er- 
rors to the real estimates of the cluster optical depths. We adopt 
a log-normal model for the errors on the T500, j estimate for the 
7th cluster: 

T500, j = T500, ; exp {ej), (25) 

with ej being a normally distributed variable of zero mean. The 
symbol fsoo, j denotes the Monte Carlo estimate of the yth clus- 
ter's optical depth obtained from the real estimate T500, y. For 
each of the 100 Monte Carlo simulations, we simulate values 
of T for each cluster and then repeat the full analysis, setting 
constraints on the kSZ monopole, kSZ-induced rms excess, and 
kSZ dipole, as outlined in previous sections. 

These analyses show that our constraints on the kSZ 
monopole, velocity rms, and dipole uncertainty change by less 



than 10% when considering errors in the optical depth estima- 
tion of the order 20-40%. Typical changes are at the level of a 
few percent (see Table 3); since we are constraining ensemble 
quantities obtained from sub-samples of the cluster catalogue, 
errors tend to average out if they are independent from cluster to 
cluster (as we expect them to be). Provided that the relative un- 
certainty in cluster luminosities is about 40 % for the L-M scal- 
ing relation (Pratt et al. 2009), and combining it with the approx- 
imate scalings L oc M'^^^ and M cc R^, one deduces that the rela- 
tive uncertainty in R should be at the level of about 10 %. For the 
spherical estimates of Tsph, 500 this translates into a roughly 30 % 
uncertainty, decreasing to about 20 % for the cylindrical optical 
depth estimate (Tcyi cc R^). From the results of our Monte Carlo 
approach above, we conclude that errors in the optical depth es- 
timates should not significantly bias our kSZ constraints. 

6. Conclusions 

The MCXC cluster sample has been used to search for signa- 
tures of peculiar velocities in the Planck CMB data. For this 
purpose, two different filters were applied: aperture photome- 
try; and the unbiased Multi-frequency Matched Filter The for- 
mer is a simple, quick and robust tool for providing estimates 
of kSZ-induced temperature anisotropics, and although it de- 
tects the tSZ-induced monopole and rms excess at the cluster 
positions, it fails to detect the kSZ eff'ect, setting a constraint on 
the kSZ-induced radial velocity rms at the level of 1200kms ' 
(95 % confidence level) for a massive and distant MCXC sub- 
sample of 1000 sources. By effectively removing the tSZ signal, 
matched filters are able to place stronger constraints, reaching 
the level of 800 km s"' (95 % C.L.) for a subsample of 100 mas- 
sive clusters. All these values, however, lie a factor of 3-5 above 
ACDM expectations for clusters of typical mass 2 x lO'"* M© at 
z^0.15. 

Both methods also provide measurements of the clusters' av- 
erage velocity that are compatible with zero: these are at the 
level of 120-160 km S-' (95 % confidence level) for the uMMF 
and AP filters. The fact that this constraint applies to a cluster 
sample whose mean redshift is z ~ 0.18 provides very strong 
evidence that the CMB is mostly at rest with respect to those ob- 
servers (as opposed to the relative motion of our local CMB to 
those sources, which is of the order of cz ~ 54 000km s"'). By 
itself, this measurement constitutes an unprecedented and valu- 
able confirmation of a prediction of the standard cosmological 
scenario, and has strong implications in discussions of the ho- 
mogeneity of our Universe. 

In this context, the large number and redshift distribution of 
the Planck cluster kSZ measurements are ideal for constraining 
void models, which attempt to explain the apparent acceleration 
without dark energy or modified gravity. Indeed, void models 
which fit the Union2 SN data are ruled out at extremely high con- 
fidence. In principle it may be possible to cancel the kSZ effect 
generated at cluster positions in these models with a large (order 
unity) isocurvature mode at last scattering ( Yoo et al. 2010), but 
this would almost certainly entail substantial fine tuning of the 
isocurvature mode. Therefore the Planck kSZ data strongly sup- 
port the conclusions of previous studies which found that void 
models generically predict very low Hq (e.g., Zibin et al. 2008; 
Bull et al. 2012) and too large kSZ power on small scales (Zhang 
&Stebbins 2011). 

Planck's constraints on the amplitude of the local bulk 
flow provide an independent view of a long ongoing debate. 
Unfortunately, our results are not sensitive to the local volumes 
where many claims for bulk flows have been raised: the limit 
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of 390 km s ' within spheres of 350/2 ' Mpc at 95 % C.L. does 
not permit us to confront claims at the level of 400-700 km s"' 
within radii of 50-120 /r' Mpc, (Hudson et al. 2004; Watkins 
et al. 2009; Feldman et al. 2010). However, on larger scales, 
Planck is able to set strict constraints on the amplitude of bulk 
flows (below 254 km s~' at 95 % C.L. for a radius of 2 /i ' Gpc), 
in clear contradiction with some previous claims (Kashlinsky 
et al. 2008, 2010; Abate & Feldman 2012). 

The linear continuity equation states that peculiar velocity 
surveys are sensitive to fluctuations in the distribution of mat- 
ter and energy on scales larger than density or galaxy surveys. 
The fact that Planck is able to set such strong constraints on 
peculiar velocities in a cluster population at (z) ~ 0.18 (and ex- 
tending out to z ~ 1), translates into correspondingly strong con- 
straints on the amplitude of primordial fluctuations at Gpc scales. 
If the Universe were inhomogeneous on scales larger than the 
size of the volume containing our cluster catalogue, these clus- 
ters would show a significant dipolar pattern in their kSZ veloc- 
ities. We conclude that Planck constraints on peculiar velocities 
are compatible with ACDM expectations, and constitute an un- 
precedented piece of evidence for the local homogeneity of the 
Universe in the super-Gpc regime. 
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Table 2. Estimated dipole coefficients (columns and velocity magnitude (column 5) using the all-sky method. The values in 
parentheses are determinations using the more restrictive mask which includes 1321 clusters. The "HFI" row reports the results 
obtained for the actual data; as discussed in the text, this estimate is significantly contaminated by Galactic foregrounds and tSZ. 
"PSM diffuse" reports the contribution from the diffuse Galactic component found in the PSM simulation. In the rows corresponding 
to thermal SZ (tSZ), instrumental noise, and CMB, the first three columns report the mean and 68 % confidence region error bars, 
while the last column indicates the 95 % upper limit. The bottom part of the table refers to results found using "shifted positions" 
for each cluster These are randomly selected between 30' and 1° from each of the MCXC clusters outside the mask region. The 
notation is as in the noise simulations, but it is relative to the distribution found for different choices of shifted positions. The velocity 
magnitude (column 5) represents the mean and 68 % error in the distribution. These rows indicate the size of the apparent dipole 
that one could find using this method, even without a cosmological dipole existing. 
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